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Abstract 

This article is motivated by the observation, that calculations of the 
Unruh effect based on idealized particle detectors are usually made in a 
way that involves integrations along the entire detector trajectory up to 
the infinitely remote future. We derive an expression which allows time- 
dependence of the detector response in the case of a non-stationary 
trajectory and conforms more explicitely to the principle of causality, 
namely that the response at a given instant of time depends only on the 
detectors past movements. On trying to reproduce the thermal Unruh 
spectrum we are led to an unphysical result, which we trace down 
to the use of the standard regularization t — > t — ie of the correlation 
function. By consistently employing a rigid detector of finite extension, 
wc are led to a different regularization which works fine with our causal 
response function. 

1 Introduction 

In 1976, W. G. Unruh d_ discovered that an uniformly accelerated observer 
moving through empty Minkowski space sees a thermal radiation with a 
temperature directly proportional to the observers proper acceleration a: 

ft 

a . 



2-Kck 



This surprising and simple result, which blurs the traditional distinction 
between "emptiness" on the one hand and matter filled space on the other, 
has given rise to numerous investigations. 
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From the many ways of treating the Unruh effect, we concentrate on 
the approach originally used by Unruh and refined afterwards by DeWitt 
[2] , which we consider the most "down to earth" , namely the concept of an 
idealized particle detector: A pointlike quantum mechanical system with dif- 
ferent internal energy states, which is coupled to a scalar field in its vacuum 
state and is moving through Minkowski space along a given trajectory x(r). 
A transition of the system from its ground state to an excited state will be 
interpreted as the detection of a particle of the corresponding energy. 

The reaction of the detector to the motion x(r) is formulated via the 
so-called response function F(lo), which gives in essence the probability of 
finding the detector in an excited state of energy u) above its ground state. 
The usual expression for this response function can be found, for example, 
in0: 

00 00 

F(u) = J dr J dr' (0|</»(x(r))^(x(r / ))|0) . (1) 

—00 —00 

It consists of a twofold Fourier transform of the correlation function of the 
field 4>, evaluated along the detector trajectory x(r). It may be considered a 
bit unfavourable, that this form of the response function does express neither 
time dependence nor causality, by which we mean the following: It has to be 
expected, that a particle detector moving through Minkowski space along a 
trajectory which is in general non stationary, will show different reactions at 
different instants of time. It is therefore desirable to have an expression for 
its reactions which allows an explicit time dependence, and, what is more, a 
time dependence of the kind, that the reaction at a given instant is influenced 
only by the part of the trajectory lying in the past of this instant. In 
there is obviously no room for such a time-dependent and causal behaviour, 
because r and t' are integrated over and F(u>) can, by construction, never 
be time-dependent. 

In section 2 we therefore consider the momentary transition rate F t (lu) 
instead of F(uj). We derive an expression for F T {u) which depends in a 
manifestly causal way on the proper time r and which can be applied to 
non-stationary trajectories. Unfortunately, the evaluation of this response 
function leads to an unphysical result already in the case of hyperbolic move- 
ment. 

In section 3 we briefly review parts of an article by S. Takagi where 
he uses a correlation function different from the standard one and tries to 
justify this by considering the pointlike Unruh detector as the limit of a 
rigid and finitely extended detector. With Takagis' correlation function, our 
causal response function leads to the expected result in the Unruh case. 

But because we consider Takagis' argumentation not entirely correct, we 
propose in section 4 our own correlation function, which we consider well 
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founded and appropriate for the problem of the moving detector. This leads 
on causal calculation to the Unruh effect without any difficulties. 

In section 5 we show an application to a non-stationary trajectory which 
smoothly starts from rest and goes over to an uniform acceleration. 

In the final section 6 we briefly state our results and formulate a critical 
remark on an often heared argument which relates the thermality of the 
Unruh radiation to the existence of a so-called acceleration horizon. 

We use the Minkowski metric with "spacelike" signature (—,+,+,+) 
and denote three-vectors by boldface, four-vectors by sans-serif letters. H 
and c are set equal to 1. 

2 Causal transition rate 

We employ as detector an idealized "atom" , i.e. a quantum mechanical sys- 
tem with a two-dimensional state space. The Hamiltonian Hp, of the detector 
atom is of diagonal form with respect to the basis |0) (ground state) and 
|1) (excited state). The energy of the ground state is 0, the energy of the 
excited state is denoted by u. The detector is moving along the trajectory 
t = t(r), x = x(t) through a massless scalar field cf)(t,x) to which it is cou- 
pled by the interaction /x(r)0(r). /x(r) is the detectors monopole moment, 
4>{t) := 0(t(r),x(r)) the scalar field at the detectors position at time r. The 
complete Hamiltonian is H = Hp, + Hp + (i(f>, where the explicit forms of 
Hp,, Hp and fi can be left unspecified. 

We now ask the following question: Suppose an energy measurement at 
time ro has given the value for the detector atom and the vacuum state for 
the scalar field. Then, what is the probability of finding the excited state of 
the detector, i.e. the energy value uj, in a later measurement at time t > tq 
if the detector has in the meantime been moving along a given trajectory? 
We do the calculation in the interaction picture, in which both, states and 
operators, are time-dependent. The field operator has the form 

*(t,x) = ^3 jSlk (^k) 6 " 1 ^^"^ W(kK Kk) <- kx) ) (2) 

where the creation- and annihilation-operators obey the covariant commu- 
tation relations 

[a(k),at(k')] = (27r) 3 2cj(k)5(k - k') , 
and uj(k) = |k|. The time evolution of the states is governed by 

4Wr))^(r)^)Wr)) 
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where 



M ( r ) = e lHDT n(0)e- iHDT . 

At To the system is in the product state \iJj(tq)) = |0)|0), i.e. both the 
detector and the field are in their ground state. Then to the first order of 
perturbation theory 

T T 

|<M<Kt)>| 2 = |(1H0)|0>| 2 | dr'j dr" e-^'-")(0|^(r')|^)^|^(T")|0) 



TO T() 



is the probability to find the detector in the state |1) and the field in the 
state | (f) at time r. Summing over the unobserved field states and using 
their completeness, we arrive at 

r r 

£|<l^(r)>| 2 = |(1H0)|0>| 2 / dr' I dr"e--^'-")(O|0(r')^»(r")|O) , 

\<p) 



TO T() 



which is the desired probability of finding the detector in its excited state 
at time r. It consists of a constant, detector-specific coefficient and an inte- 
gral which depends only on the field (j). We are interested only in the inte- 
gral which we denote by F t (uj), and we introduce the abbreviation W(t,t') 
( "Wightman-function" ) for the correlation function (0|</>(t)</>(V)|0): 

T T 

F t {uj) = J dr' J dr" e-^'-^W^r') . 

TO T 

The integration extends over a square in the r', r"-plane. We can write this in 
a more convenient form by introducing new coordinates u = r' , s = t' — t" 
in the lower triangle t" < r, and u = r" , s = r" — t' in the upper triangle 
t" > t'. This leads to 

T U — To 

F T {u) = j du j ds (e~ iuJS W(u,u- s) + e +iuJS W (u - s , u)) , 

T 

which is the same as 

r u— to 

F t (uj) = 2 J du J ds Re (e~^ s W(u, u - s)) 

TO 

because of W(t',t) = W*(t, t'). We can read off the following structure: 
The probability of finding the detector at time r in a state, whose energy is 
different by an amount u from the energy of the ground state, is given as a 
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sum of contributions from every time u between To and r. The contribution 
of a certain time u is in turn given by the correlations of <fi(u) with (j) at 
every past time since tq. Therefore, the time derivative of F t (uj), i.e. the 
transition rate F T (co), is not determined only by the state of motion at time 
r but by the whole course of the trajectory between tq and r: 



T— TO 

F t (lo) = 2 J ds Re (e~ iujs W(t, t - a)) 



From now on we will consider only the transition rate F t (uj) which is pro- 
portional to the number of "clicks" per second in a detector consisting of a 
swarm of identical detector atoms, and which is physically more accessible 
than the probability F T (uS). (In the following we will call F T {uo) "response- 
function", which normally rather is the name for the probability F t {uj).) Fur- 
thermore we get rid of the arbitrariliy chosen tq by doing the limit tq —* — oo: 

oo 

F t (uj) = 2 J ds Re (e~ iuJS W(t, t - s)) . (3) 
o 

If W is invariant under time translation, i.e. W{t\ + At, T2+ At) = W(t±, T2) 
then it is an easy matter to show that 

00 00 

2 j ds Re (e' iujs W(r, T-a))= j ds e~ iuJS W(s, 0) . (4) 

-00 

That means, if W is invariant under time translation, F t (uj) is independent 
of r and can formally be written as an integral over the whole (i.e. past 
and future) trajectory. This does not lead to any problem in the case of 
a stationary trajectory, because this is a "global" object anyway, for an 
arbitrarily small piece of such a trajectory contains all the information about 
its whole, in particular its future, course. But in the general case of a non- 
stationary trajectory, one must use (jSJ) which is therefore the starting point 
of our further considerations. 

Up to now, our reasoning was completely general and we did not make 
any use of the fact that <ft is a scalar field of the form ©. In the infinite- 
dimensional case of the scalar field, it has to be taken into account that W 
is not a function but a distribution. This is clear from 

1 r d 3 k d 3 k' , 1 1 \ 

<0«x)«x<)|0) _ _ ( „|y__ a(k)a t (kV .(--)|0) 

e ik(x-x') 



(2tt) 3 J 2cu(k) 
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because the integrand is oscillatory and the integral does not converge in the 
classical sense. The standard way to handle this is to make the replacement 
t — > t — is, or, equivalent, to insert an exponential factor e~' k ' £ into the 
Fourier representation, which regularizes the high-frequency behavior and 
take the limit e — * at the end of the calculation, i.e. 

WWMIO) = ^ J ^ e*<*-"H"l< (5) 

and 

oo 

F t (uj) = 2 \iin J ds Re (e~ iuJS W £ (T,T-s)) . (6) 
o 

We will investigate later how this regularization can be done in a more 
physically motivated way. Accepting the e~' k ' £ -factor for the moment, one 
arrives after a short calculation at 

w(xwxoio )= (t _ t ,_-y: 2 (x _ x , )2 ■ (t> 

This is the form of the correlation function which is used in the standard 
literature, see e.g. J3j, (3.59). If the trajectory of a detector at rest, t(r) = 
r, x(t) = is inserted, one gets 

oo oo 

F t {uj) = — r— 9 lim / ds Re — j = lim / ds 



2tt 2 e^o J (s - is) 2 4tt 2 e—o J (s - is) 2 

-oo 

One can evaluate this integral by contour-integration along the real axis 
closed by an infinite semi-circle in the upper or lower half-plane, depending 
on the sign of u. The resulting value is 2iruje £UJ 0(— oS) which leads to 

F T (u) = -— lj@(-u) . 

This result - independent of r as expected - describes the fact that a detector 
atom at rest does not go spontaneously into its excited state (F t (uj) = 
Voj > 0), but of course it decays spontaneously into its ground state due 
to the coupling to the scalar field. 

Now the interesting thing is to see what happens if the detector is uni- 
formly accelerated, i.e. 1 

t(r) = sinh(r), x(t) = cosh(r), y(r) = z(t) = . 



1 We set the proper acceleration equal to 1 because the transition-rate for arbitrary 
acceleration a > can be obtained by a simple scale-transformation F t (lo) = aF^-~ r> ( — )■ 
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Because this is an orbit of the Killing vector xdt+td x it is stationary and one 
expects a time-independent transition rate. Therefore it is a bit surprising 
to see that insertion of this trajectory into ® leads to a W e (r, r — s) which 
is not time-independent, i.e. not independent of r: 

WAr.r-s)- - 1 ' 4 " 2 



4sinh (I) — e 2 — 2ie (sinh(r) — sinh(r — s)) 



But this in itself is not to alarming, because time independence is demanded 
only for the transition rate, which contains a limit e — > 0, and Wq(t,t — s) is 
clearly independent of r. To settle the case, let us do a numerical calculation 
of F T {u)): The response function © is of the form 

oo 

F t {oj) = 2 lim / ds (a £ (r, r — s) cos us + b e (r, r — s) sin us) , 
e-^OJ 


where a £ := ReW £ and b £ := lmW £ , i.e. 

a £ (r,T-s) = 

b e (r,T-s) = 



4sinh 2 (f) -e 2 



47r2 (4sinh 2 (f ) - e 2 f + 4e 2 (sinhr - sinh(r - s)f 
1 2e(sinhr — sinh(r — s)) 



47r2 (4sinh 2 (f ) - e 2 f + 4e 2 (sinhr - sinh(r - s)) 2 

In order to calculate the transition rate for a given time r, we thus have to 
evaluate the cosine- or sine-transform of a e (r, r — s) and 6 e (r, r — s) with 
respect to s at fixed r. The numerical results for seven different times r = 
0, ±1, ±2, ±4 are shown in figure^ The lower curve corresponds to r = —4, 
higher curves to increasingly later times. Thus our suspicion is confirmed 
that the r-dependence of the correlation function leads to a r-dependent 
transition rate, which is unacceptable in the case at hand. Furthermore, the 
results are problematic for r < because F t (oj) assumes negative values in 
this cases. 

At this point, the first question which comes to mind is to ask the relia- 
bility of the numerical scheme. Fortunately, an exact calculation is possible 
at least in the two special cases of r = or u = 0. In the case r = the 
transition rate can be evaluated by the method of contour-integration. This 
results in 

P ^ = ^^—l > 

which is indeed the usual Unruh spectrum. This is exactly in agreement with 
the numerical result for r = 0. In the other case, u = 0, the integral can be 
evaluated exactly with the help of Maple, which leads to 

l + e 2r (2r- 1) 



F T (0) 



2ir 2 (e 2T -1) 
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Figure 1: Numerical results for F T {oo) 



In this way we get on the one hand a rigorous proof that F t {uj) is indeed 
r-dependent if the calculation is based on ((7J). On the other hand this exact 
result permits a further test of the numerical calculation. For example, the 
numerics tells us that -Fi(O) = 0.0104..., and indeed 

For any other checked F T (0) there is an equally good agreement between the 
numerical prediction and the analytic result. This gives us confidence in the 
correctness of the numerical calculation in the general case. 

We are thus confronted with the following situation: The application of 
our causal formulation of the detector response to the case of uniformly ac- 
celerated motion lead to a unphysical, time-dependent result. But we see no 
reason to believe that it is impossible to derive the Unruh effect consistently 
in an explicitly causal way, or, more generally, to formulate a causal response- 
function which can be applied to non-stationary trajectories. Therefore, we 
will seek the error not in the use of our causal response function but 
rather suspect that the correlation function (JJJ) and the e~' k ' E -regularization 
of the integral representation (jSJ) is not suitable for the problem of the mov- 
ing detector. In the following sections we will consider a finite, rigid detector 
and show that one is led in a natural way to a correlation function which is 
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different from 2 . We use the article (3] of S. Takagi as a model, even if we 
come to a different result. The new correlation function leads upon causal 
calculation to the desired Unruh result without any difficulty, and it seems 
to be superior to (J7J) in other ways. 

3 Correlation function a la Takagi 

In [3] Takagi criticizes the ad-hoc regularization of the correlation function 
by the introduction of the factor e~' k ' £ into its Fourier representation or 
by the equivalent prescription t — > t — ie. Rather he would like to understand 
the origin of the regulator as a consequence of the finite extension of the 
detector. Therefore he uses instead of the ultralocal 4>(t) = 0(x(r)) the 
"smeared" version 

<f>(r) = J d 3 e/(O0(x(r,O) , 

where /(£) is a weight function which is normalized, symmetric under ro- 
tations and concentrated around the origin. By choosing /(£) = <5(£) one 
recovers the pointlike detector; we will consider a family of finite detectors 
/ e (£) instead with lim e ^o /e(£) = <^(£)- By taking the limit only at the end 
of the calculation, we thus model an "infinitesimal" instead of a pointlike 
detector. 

Of special importance is the function x(r, £), so we have to say a little 
bit about it. It describes the transformation from the usual Minkowski coor- 
dinates to the Fermi- Walker coordinates r and £:=(£, f], C) associated with 
the trajectory 3 x = x(r): 

*(t,0 = *(t) + £f(t) 

x(r,0 =1^1. (8) 

The geometrical meaning of these coordinates is the following: To every 
point x(t) = (t(r), x(r), 0, 0) on the trajectory belongs its simultaneity, i.e. 
the hyperplane orthogonal to the four- velocity u(r) = (i(r), x(t), 0, 0) in this 
point. This three-dimensional space consists of all events which are simulta- 
neous to x(t), where simultanity is judged from the co-moving inertial frame. 
If we attach orthonormal basis vectors e^(r) = (x(t), t(r), 0, 0) , e r? (r) = 
(0, 0, 1, 0) , (r) = (0, 0, 0, 1) to every such hyperplane, we can characterize 
every event x (in a neighborhood of the trajectory) by (r, £,77, £): First, we 

2 In HH the Unruh spectrum is derived from J7J only by the use of a not entirely con- 
vincing calculation in which "a positive function of r, r' is absorbed into e". 

3 For the moment, we restrict ourselves to trajectories which run along the a:-axis, i.e. 
y = z = 
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fix r by the condition that the given spacetime point lies in the simultan- 
ity of the trajectory point x(r). Then x — x(r) is a linear combination of 
e^(r), e^(r), e^(r), and we can read off (£,r/, () as the corresponding coeffi- 
cients. That is 4 

x(t) + £e 5 (r) + t]& v {t) + Cec(-r) = te t + xe x + ye y + ze z , 

which is nothing else but (|HJ)- Expressed by the coordinates (r, £), the Minkowski 
metric assumes the form 

ds 2 = - (1 + 2(ix - xi)£ + (x 2 - dr 2 + d£ 2 . 

Obviously, every simultanity r = const, carries the same cartesian three- 
metric. Because <P(t) originates from an integration over the r-simultanity, 
the r- independent shape function ensures that this construction cor- 
responds indeed to a rigid detector - where rigidity means that its three- 
geometry as seen from its own momentary rest system is unchanged in the 
course of proper time. In the case of hyperbolic motion we have for example 

x(r,£) = 

and the metric 

ds 2 = -(l + £) 2 dr 2 + d£ 2 , 

i.e. the well known transformation to Rindler coordinates, cf. 0. 
Now by linearizing x(r, £) around £ = 

x(t,0 «x(t) + , „ =: x(t) + ^(r) 




Takagi finds 

0(T) = (2^)3 
where 



/ (°( k ) e " ikx(T) />;r) + at(k)e ikx «f(k;r)) , 

/(k;r) := | d 3 £ /(£) e - k «W 

converges to for |k| — > oo, which makes 0(t) a well-defined operator. He 
now chooses for / (£) the family of functions 

1 ^3 



7T 3 (e + e 2 )(r/ 2 + e 2 )(C 2 + e 2 ) 



4 r is used with two different meanings: As proper time of the trajectory, and as space- 
time coordinate. This should not lead to any confusion. 
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which converges to the three-dimensional (5-function in the limit e — > 0. Then 
/(k; r) can be calculated explicitly: 

/ £ (k;r) = e~ £|kX « (r)l ~ £|kX " (r)h£|kX « (r)! . (9) 

Takagi now claims that this complicated expression can be effectively re- 
placed by the simpler 

/ £ (k;r) = e -l k l £ / 2 , (10) 

which does not depend on r any more. Accepting this for the moment, one 
arrives after a few lines at 

(oWtMtOIo) = (2^)3 / ^ eik(x ~ x ' Hk|£ • ( u ) 

Obviously, the finite extension e > of the detector manifests itself in a 
regularization of the integral, as desired. Takagi now writes this in four- 
dimensional form 

(O|0(r)0(r')|O) = JL_ Jd'k e(k°)5(k 2 y^-^ k °c , (12) 

and claims this form to be manifestly Lorentz invariant. He now uses this 
Lorentz invariance to derive 

(O|0(r)0(r')|O) = . (13) 



sgn(i - t') ^/-(x-x') 2 - ie) ' 



By inserting the hyperbolic trajectory and using contour-integration, it fol- 
lows indeed that 

• 1 to 

Even if this is the desired result, we have some reservations concern- 
ing the derivation. First, the claim that the integral ((12)1 is Lorentz invari- 
ant holds only if e = 0. But in our case e > the factor e _fc ° e destroys 
this invariance, which makes the evaluation of the integral obsolete. Indeed, 
straightforward evaluation of 1)1 ljl leads exactly to the standard form of 
the correlation function and not to 1)13)1 . Therefore, if Takagi's derivation 
were correct up to 1)11)1 this would show that taking into account the finite 
extension of the detector does not solve the problem which we encountered 
in the last section. We now suppose, that, contrary to Takagi's assertion, it 
is not allowed to replace © by 1)10)1 . Indeed this replacement eliminates the 
dependence on r and on the trajectory, and it is not clear that this leaves 
the final result unchanged. Therefore we will try to work out Takagi's idea 
more carefully in the next section. 
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4 Correlation function of an infinitesimal rigid de- 
tector 



We now use the monopole interaction with a "smeared" field-operator 

<j>{r) = J d 3 e/(0^(x(r,0) (14) 
instead of <j)[r) = </>(x(r)). The correlation function thus assumes the form: 

(O|#(t)0(t')|O) = J^{0\J d'{ /«) J d'Y /«') 

Using the commutation relations of the a and a^, a|0) = and (0|0) = 1 we 
can do the k'-integration and find 

««t)#(t')io) = _i_y_^_ 9(k;T)9 . (k;Tl) 

where we have introduced g(k;r) := J d 3 £ /(£)e* kx ( r, £). In order to find the 
correlation function of a detector with shape /(£) and trajectory x(r) we 
thus have to determine the transformation to Fermi coordinates x(r, £), from 
that the function g(k;r), and finally (0\4>(t)4>(t')\0). This amounts to six 
integrations, so one would expect difficult computations for even the simplest 
trajectories. In fact we will see, that by choosing a suitable shape / (£) of 
the detector, one arrives at an explicit form of the correlation- function for 
any arbitrary trajectory. 

Let us first consider the simple example of a detector which is not moving 
at all, and whose shape is given by a function from the family 

f^ = hw^? ' £>0 • (15) 

The / e (£) are invariant under rotations, satisfy f d 3 £ f £ (£) = 1 Ve > and 
have the scaling property f £ (£) = £~ 3 /i (£/£)• This means, they approximate 
the three-dimensional (5-function. The parameter e has the physical dimen- 
sion of a length and indicates the extension of the detector: / e (£) is negligible 
outside £ 2 < e 2 . By taking the limit e — > at the end of the computation 
one can model an infinitesimal detector. The transformation to Fermi coor- 
dinates is trivial in the case of a detector at rest: i(r, £) = r, x(r, = 
This leads to 

g £ (k; t)= J d 3 e fs(ty kx(T ^ = e-^ e -sM 
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From that: 



1 f c\ 3 k 

(0|«Kt)<Kt')|0> = ^3 y ^fe(k;r) S :(k;T') 

1 f e -*|k|(T— r'-te) 



d 3 fc 



(2tt) 3 7 2|k| 

-l/4vr 2 
(r — r' — ie) 2 



(16) 



(We have replaced 2e by e.) Using the same contour-integration as before in 
section 2 it follows 

F T (u) = -—uQ(-u) . 

So in this simple case we get the same result as before in section 2 where we 
inserted the trajectory into the Wightman function Q. But we now see that 
the high-frequency regulator e~ e ' k ', which was put in by hand in section 2, 
results from the finite extension of the rigid detector and thus acquires a 
physically meaningful interpretation. 

Let us now consider an arbitrary trajectory x = x(r). The transformation 
to Fermi coordinates is 

x(r, £) = x(r) + ^(t) + r/e^(r) + Ce c (r) , 

where e^(r), e r? (r), e^(r) together with u(r) = x(r) form an orthonormal 
basis which is Fermi- Walker transported along x(r). We will use the same 
detector-shape (fTT))) as in the case of the detector at rest. Let us start with 
the calculation of g e (k;r): 

5e (k;r) = J d 3 ^ / £ (£)e ikx(T ' = e* kx{r) J d 3 £ /^e^M^+^M+^cM) 

We now introduce the abbreviation k := (ke^(r), ke ?? (r), ke^(r)). It is easy 
to do the integration by using spherical coordinates in ^-space: 



That means 

5e (k;r) = e ^xM e -£ V / ( ke «( T )) 2 +( ke ^ T )) 2 +( ke cM) 2 

Because k is a lightlike four- vector, the sum of squares under the square root 
is a complete square: By multiplying the identity 

k = -(ku)u + (ke 5 )e f + (ke ?? )e r? + (ke^)e^ 
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with k we get 

(ke € (r)) 2 + (ke„(r)) 2 + (ke c (r)) 2 = (ku) 2 = (-|k|t + kx) 2 . 

ku is the product of a lightlike and a timelike vector and therefore negative, 
so that 

5£ (k;r) = e ikx(r) e +£ku(T) . 

Obviously, e eku ^ and not simply e~ £ ' k ', is the physically motivated regular- 
isating factor resulting from the finite extension of the detector. Using this, 
we now calculate 

Jk(x-x') e +ek(u+u') 



J TjX fc(k;T) 5 *(k;r') = J ^- e l 



+ 



From this we finally obtain 



., . , . .. .... e -i\m-t'-Hi+i')) 

^3^, ik(x— x— is(x+x )) 

|k| 

47T 

-(t - t> - ie(i + i')) 2 + ( x - x' - ie(x + x')) 2 

47T 



(x — x' — ie(x + x')) 



2 



W(r)^)|o)= r^TTTTT^ • (17) 

(x(r) - x(r') - z£(x(t) + x(r'))J 
This result has to be compared with 

(Q|^(x(r))0(x(rO)|O) = , (18) 

(t(r) - t(r') - te) - (x(r) - x(t')) 

which is obtained from the usual Wightman function (J7J) by inserting the 
trajectory x = x(r). It is remarkable, that (|T7|) . in contrast to cannot 
be written as a function on Minkowski space because it depends on the 
four-velocity x. It can easily be seen, that (|18jl results from (|17j) by inserting 
x(t) = (1,0,0,0), which, of course, is inconsistent because the detector is 
not at rest but moving along the trajectory x = x(r). As a first test of our 
new correlation function, we insert all the possible stationary trajectories in 
Minkowski space, which are classified in into six groups. It turns out that 
in all six cases the correlation function depends only on r — r' which has to 
be expected for a stationary movement. In contrast, inserting into (|18|) leads 
to a result which is invariant under time translation only in the cases of the 
detector at rest and in circular movement. We see this as further evidence 
that lfT%|) is not suitable for the calculation of the detector response and 
that (|17[) should be used instead. Furthermore, (J17|) is manifestly invariant 
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under Lorentz transformation of the trajectory, which we consider as another 
advantage over Q18|) . For example, inserting t(r) = jt, x(t) = 7vr with an 
arbitrary v satisfying v 2 < 1 leads to the same correlation function IJ16JI as 
in the case of a detector at rest. 

Concluding: Our proposal for calculating the transition rate at time r 
for a detector moving on x(r) is 



F t (lo) = 21im^ds Re (e- iuJS (0\<P(t)<P(t - s)|0» , (19) 
o 

where 

WW(tO|0) = — ■ (20) 

(x(r) - x(r') - ze(x(r) + x(r'))j 

Thus, in order to calculate the time-dependent response of a detector in a 
manifestly causal way, i.e. by using only information on its movements in 
the past, we use a different integration formula and a different correlation 
function. 

It is now easy to show that in the case of an uniformly accelerated 
trajectory the correlation function is 

<oMW)|o> = - ~ 1/167r2 - — 2 , 

(sinh (^f) - ie cosh {^f)) 

which clearly is invariant under time translations, as has to be expected from 
the stationarity of the trajectory. We can thus write the detector response 
with an integration over the whole real axis, according to iflj): 

oo 

F t (lu) = — — Q lim / ds k 

167T 2 ^0_7 (sinh(f) -^cosh(f)) 2 

and evaluate the integral with the help of contour-integration: We integrate 
along the real axis in positive direction and along the line Ims = 2tt in nega- 
tive direction. The only singularity within this contour lies at s = 2i arctan e. 
A short calculation shows that because of the periodicity of cosh and sinh, 
the value I £ (w) of the integral is connected to the value of the residue by 

( e -iu)s \ 
2 ' s = 2i arctan e = I £ (lo) - e 27TUJ I e (uj) . 
(sinh (f) -ie cosh (f)) / 

The value of the residue is found to be — Aiuje 2u) arctan e /l + e 2 which leads to 

8itlo e 2ujaictan£ 
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and finally 

■ 1 ijj 

This is of course the well-known thermal Unruh spectrum, which we derived 
here in a different way as usual, namely in a time-dependent, causal manner 
and by using a physically founded regularization of the correlation function. 
This approach may be considered more satisfying than the usual one. 

5 Application to a non- stationary example 

Let us consider the non-stationary trajectory 



t{r) = r-ln2 + yi + ie 2 ^-ln M+^l + I e 2rJ 

x(t) = ~e r , y(r) = , z(t) = . (21) 

This movement interpolates smoothly between rest at x = and uniform 
acceleration along x 2 — t 2 = 1. Actually, it is easy to show that t(r) — ► sinhr 
and x(t) — > coshr if r — > +oo. For r — ► — oo we have instead i(r) — * 
t + 1 — 2 In 2 and x(t) — ► 0. The square of the four-acceleration 

h 2 = _p +i 2_ 1 



1 + 4e~ 2T 



shows that the proper acceleration a = vb 2 is indeed smoothly increas- 
ing from to 1. The radiation spectrum is now time-dependent, of course. 
The result of a numerical calculation for three times r = —4, 1,6 is shown 
in figure [21 At time r = — 4 the movement has not really yet begun, so 
one gets the spectrum of the detector at rest. At time r = 6 the move- 
ment ran a fairly long time with approximately constant acceleration a = 1 , 
so the thermal spectrum results, which does not change any more at later 
times. It is remarkable, that at the intermediate time r = 1, i.e. a = 0.65, 
the transition-rate decays slower for high to than in the case of the ther- 
mal spectrum at late times. This means, that the increasing of acceleration 
takes place smoothly but not adiabatically: The radiation spectrum is of a 
non-thermal nature, and not simply thermal with only a time- varying tem- 
perature. This shows that not only acceleration contributes to the radiation 
effect, but higher time derivatives as well. 

6 Conclusion 

Starting with the intention to derive an explicitly causal formulation for the 
transition rate of a non-inertial Unruh-DeWitt detector, we were led to the 
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Figure 2: Non-stationary Unruh-like spectrum 



question of which regularization of the scalar field's correlation function is 
appropriate to this problem. The usual prescription t — > t — ie or the equiva- 
lent insertion of the convergence factor e _£ ' k ' into the Fourier representation 
was shown to be inappropriate by a numerical calculation, which lead to a 
time-dependent unphysical result in the case of uniformly accelerated mo- 
tion. Following Takagi we arrived at the correct regularization by considering 
the pointlike detector as the limit of a finite, rigid detector - where rigidity 
is defined with respect to its own momentary rest system. The result is 

oo 

F T (u) = 2 lim J As Re (e~ iuJS (O\0(t)0(t - s)|0)) , 

l/4vr2 

(x(r)-x(r')-ie(x(r)+x(r'))) 2 

We consider this as the physically well-founded form of the response function 
of an infinitesimal detector. 

Finally, a short criticism of an often heared argument regarding the 
Unruh effect may be formulated. The fact that Unruh radiation is exactly 
thermal is astonishing and demands for some deeper explanation. A well- 
known argument is based on the existence of an acceleration horizon which 
seperates the detector from some degrees of freedom of the scalar field. 



o 

W(t)<Kt')|0) = 
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Indeed, the eternal uniformly accelerating detector or observer is causally 
separated from the "left half" of Minkowski space by the horizon t = x. 

Therefore, as the argument goes, the observer describes the global vac- 
uum state |0) of the scalar field as a density matrix which results from 
tracing out the degrees of freedom behind the horizon. This density matrix 
turns out to be thermal, with a temperature of T = a/27r, cf. [Jj. In this way 
the ocurrence of the thermal spectrum can be understood as a consequence 
of the existence of an acceleration horizon and the associated "hiding" of 
degrees of freedom. This may look convincing, but is not without problem 
if considered from a causal viewpoint: On the one hand, the existence of an 
acceleration horizon can be stated only if the whole course of the trajectory 
is known, in particular the future behaviour up to t = +00. The existence of 
a thermal spectrum, on the other hand, is a fact which results only from the 
past movement of the detector. Then the problem is posed, that an observer, 
who percives thermal radiation here and now, cannot explain this by show- 
ing up the acceleration horizon, because this is not an observable border, 
existing somewhere in spacetime, but rather an abstract object which can 
be constructed only "posthumously" after the entire movement is known. 
At any finite time, there exist for the accelerated observer no more and no 
less unobservable degrees of freedom as for the observer at rest. 

The problem which - in my opinion - really calls for an explanation, 
is the fact, that an observer moving along trajectory (|21[). which smoothly 
goes over from rest to an arbitrarily long period of almost uniform acceler- 
ation, sees a radiation spectrum which is thermal to any arbitrary precision 
- without making any assumptions on the future movement of the detector. 
Of course, the trajectory (|21|) is defined for all t, and with it comes the accel- 
eration horizon t = x, but this does not enter the calculation. The spectrum 
measured at time r would not change at all if the detector returns to rest 
afterwards. Thus we have the following situation: A detector which is asymp- 
totically at rest for t — > ±00, is moving for an arbitrarily long (but finite) 
time with almost uniform acceleration and percives an (almost) thermal 
radiation-spectrum. Because he returns from asymptotic rest to asymptotic 
rest, there is no acceleration horizon in this case, and no hidden degree of 
freedom. The question for a simple explanation of the (arbitrarily precise) 
thermality of the observed radiation is but as important in this case as in 
the idealized case of eternal hyperbolic movement. 
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